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ABSTRACT 

Recently a second type of spicules was discovered at the solar limb with the Solar Optical Telescope 
(SOT) onboard the Japanese Hinode spacecraft. These previously unrecognized type II spicules are 
thin chromospheric jets that are shorter-lived (10-60 s) and that show much higher apparent upward 
velocities (of order 50-100 km s~^) than the classical spicules. Since they have been implicated 
in providing hot plasma to coronal loops, their formation, evolution and properties are important 
ingredients for a better understanding of the mass and energy balance of the low solar atmosphere. 
Here we report on the discovery of the disk counterparts of type II spicules using spectral imaging 
data in the Ca II 854.2 nm and Ha lines with the CRisp Imaging SpectroPolarimeter (CRISP) at the 
Swedish Solar Telescope (SST) in La Palma. We find rapid blueward excursions in the line profiles 
of both chromospheric lines that correspond to thin, jet-like features that show apparent velocities 
of order 50 km s~^. These blueward excursions seem to form a separate absorbing component with 
Doppler shifts of order 20 and 50 km s~^ for the Ca II 854.2 nm and Ha line respectively. We show 
that the appearance, lifetimes, longitudinal and transverse velocities and occurrence rate of these 
rapid blue excursions on the disk are very similar to those of the type II spicules at the limb. A 
detailed study of the spectral line profiles in these events suggests that plasma is accelerated along 
the jet, and plasma is being heated throughout the short lifetime of the event. 

Subject headings: Sun: atmosphere — Sun: chromosphere — Sun: magnetic fields — Sun: atmospheric 
motions 



1. INTRODUCTION 

Spicules are slender features protruding from the solar 
limb when observed in for example the Ca II H and Ha 
lines. Classically, spicules have been reported to show 
velocities of order 20-30 km s~^, lifetimes of order 5-10 
minutes and hei ghts up to 5,00 0-10 ,000 km above the 
photosphere (see Beckers (1968) and Sterling (2000) for 
reviews on the older literature on spicules). 

The high spatio-te mporal resolutio n of the Solar Opti- 
cal Telescope (SOT, [Tsuneta et al.|[2Q08| ) onboard Hin- 
ode dKosugi et al.||2007|) have revolutionized our view of 



spicules. 'De Pontieu et al 



(2007b| showed that spicules 
can be grouped in two categories based on Hinode/SOT 
Ca II H observations of the solar limb. Type I spicules ap- 
pear to rise up from the limb and fall back again. These 
structures show a simi lar dynamical evolution as active 
region dynam ic fibrils ([Hansteen e t al.|[2006[ De Pontieu 
et al.l 2007a )_and a subset of quiet sun mottles (|Rouppe 
van der Voort et al. 2007), so they can be interpreted as 



their ott'-limb counterparts. Numerical simulations based 
on radiative MHD indicate that the motion of type I 
spicules is driven by magneto-acoustic shocks that can 
be generated by a variety of processes, such as convec- 
tive buffeting, granular collapse, dissipation of magnetic 
energy, and leakage of waves into the chromosphere (e.g., 
Hansteen et al.|[2QQ6| [De Pontieu et al.||2007a( [Heggland 
et al.||2007| |Martinez- Sykora et al.||2QQ9 ^ 
Type II spicules show different behavior. They exhibit 
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upward motion, followed by rapid fading from the Hinode 
Ca II H passband, without a downward moving phase. 
Sometimes they accelerate while they rise. They have 
lifetimes between 10 and 100 s, apparent velocities be- 
tween 50 and 150 km s~^ and widths between 150 and 
700 km and undergo a swaying motion caused by the up- 



ward p ropagation of Alfvenic waves (De Pontieu et al 



2007c). These spicules are longest in coronal holes, with 



heights up to 10 Mm, while they appear shorter in ac- 
tive regions, where they rarely extend more than 2 Mm 
in height. 

Many questions about type II spicules remain. For ex- 
ample, it is unclear whether the apparent motions that 
are measured at the limb (which are based on the tempo- 
ral evolution of Ca II H intensity in the plane of the sky) 
are asssociated with bulk mass motion. Previous studies 
have found plentiful evidence for line of sight motions of 
order 20-30 km s~^, but little evidence for ffows of or- 
der 50-100 km s~^. Clearly, measuring Doppler shifts on 
the disk could shed light on this issue. However, while 
type II spicules are ubiquitous at the hmb, both in coro 



nal ho les and in quiet sun regions (De Pontieu et al. 



2007b), it has been unclear what their counterparts are 
on the solar disk. Finding such a disk counterpart is cru- 
cial for several other reasons. It can help reveal what 
the formation mechanism of these jets is (which is diffi- 
cult to study at the limb because of the line of sight su- 
perpositon). Type II spicules have also bee n implicated 
in pro viding the corona with hot plasma by |De Pontieu] 



et al. ((2009) who established a correlation between faint 
upffowing signals at coronal temperatures (from asym- 
metries in EUV spectral lines) and the highly dynamic 
signals in the Hinode/SOT Ca II H passband. The latter 
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signals are thought to be associated with type II spicules, 
but unequivocally establishing the disk counterpart of 
type II spicules would help reduce the uncertainties in- 
volved in establishing the importance of type II spicules 
in the mass and energy balance of the corona. 

In this paper we will focus on establishing what type II 
spicules look like on the disk. This is not a straightfor- 
ward task, given the high spatial and temporal res o lution 



required to observe these features. Wang et al.l ( 1998 



he disk 



and Chae et al. (1998) looked for lia jets on the dish 
with a spatial resolution of 1 arcsec. They found roundish 
darkenings in the blue wing of the line, with typical sizes 
of 3-5 arcsec and lifetimes of 2 minutes, without subse- 
quent redshift. These darkening are often associated with 
converging magnetic dipoles at supergranular boundaries 
in the photosphere, even though they also observe some 
jets in unipolar regions. Using cloud modeling they find 
typical upflow velocities of 30 km s~^. While some prop- 
erties fit with the observed behavior of off-limb Ca II H 
type II spicules, their sizes and lifetimes indicate that 
these are not the same phenomenon . 



More recently. 



([2008) investigated 



Langangen et al. 
so called "rapid blueshifted events (TlBEs) in on-disk 
Ca II 854.2 nm data obtained by the Interferometric 
BIdimensional Spectrometer (IBIS). RBEs are a sudden 
widening of the line profile on the blue side of the line, 
without an associated redshift. The RBEs are located 
around the network - but not directly on top of individ- 
ual network elements, show blueshifts of 15-20 km s~^, 
and have an average lifetime of 45 s. The authors inter- 
pret these RBEs as chromospheric upflow events with- 
out subsequent downflow and suggest that these might 
be the on-disk counterparts of Type II spicules. Using 
Monte Carlo simulations they show that the low observed 
blueshift can be explained by a wide range of spicule ori- 
entations combined with a lack of opacity in the upper 
chromosphere. 

Full 3D non-LTE radiative transfer calculations on a 
model atmosphere b ased on a snapshot o f radiation- 
MHD simulation by [Leen aarts et al.| ( 1200 9) show that 
darkenings in the blue wmg of the Ca II 854.2 nm line 
indeed correspond to chromospheric upfiows. 
The evidence for the connection between RBEs and 



type II spicules presented by Langangen et al. (2008) is 
tantalizing but not conclusive! With their fast evolu- 
tion and small spatial dimensions, type II spicules are 
an elusive phenomenon that are on the limit of what is 
observable with p resent day telescopes. In this study 
we improve on the Langangen et al. ( 2008| observations 
in terms of spatial and temporal resolution, and em- 
ploy a spectral sampling that is better suited to cover 
high Doppler shifts. In addition, we observed both the 
Ca II 854.2 nm and the Ha spectral lines. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were obtained with the CRisp Imag- 
ing SpectroPolarimeter (CRISP, Scharmer et al. 2008) 
installed at the Swed ish 1-m Solar Telescope (SST, 
Scharmer et al. 2003a) on La Palma (Spain). CRISP 
is a spectropolarimeter that includes a dual Fabry-Perot 
interferom eter (FPI) system similar to that described by 
2006). CRISP is equipped with 3 high-speed 
cameras that operate at a frame rate of 



Scharmer 
low-noise 

35 frames per second and an exposure time of 17 ms. 
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Fig. 1. — Details from one RBE observed in Ho;. The top panel 
shows the spectral evolution (Af-dataslice) in one spatial location 
in the RBE structure. The location is marked with a cross in the 
Ha blue- wing image (Av = —59.4 kms"-*^) shown in the bottom 
panel. A detailed line profile is shown in the middle panel. The 
profile with open diamonds is the Ho; profile averaged over the 
whole observed FOV, with the diamonds marking the sampling 
positions. The dashed line is the difference of the average and 
RBE profiles. 

The 3 cameras are synchronized by means of an optical 
chopper, 2 of these cameras are positioned behind the 
FPI after a polarizing beam splitter, and the 3rd cam- 
era is positioned before the FPI but after the CRISP 
pre-filter. The latter camera is used as anchor channel 
for image processing and is referred to as the wide-band 
channel. The image scale is O^OTl/pixel and the field of 
view (FOV) ^71x71 arcseconds. CRISP allows for fast 
wavelength tuning (^50 ms) within a spectral region and 
is ideally suited for spectroscopic imaging of the chromo- 
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Fig. 2. — Images from the dataset of 2008 June 15. Top left: blue wing of Ha line at Av = —45 km s~^; Top right: Ha line center; 
bottom left: Fe 630.2 nm blue wing (—4.8 pm) Stokes V/I^^qj^^^ magnetograrn, linear lyscaled between ±7%; bottom right: Ca H 854.2 nm 
line core. The black boxes indicate regions of interest (ROIs) shown in Figs. Is] [6]and|7| The yellow and green lines show detected paths of 
respectively Ca H 854.2 nm and Ha RBEs. The RBEs along the red and blue nnes are shown in detail as images in Figs, [s] and [7| and as 
time-slices along their trajectories in Figs. [4] and [s] Note that the Ha line core image shows significantly more fibrilar structures than the 
Ca n 854.2 nm image, indicating that the Ha line has the largest chromospheric opacity of the two. 



sphere where the dynamical evolution time can be on the 
order of a few seconds, sometimes even faster than 1 s 



van Noort & Rouppe van der Voort 2006). For Ra the 
transmission FWHM of CKiSF is 6.6 pm, and the pre- 
filter FWHM is 0.49 nm. For Ca II 854.2 nm the trans- 
mission FWHM is 11.1 pm, and the pre- filter FWHM is 
0.93 nm. 

In this paper, we analyze two data sets from 2008 June 
15, from an area at the edge of a coronal hole close to 
disk center at {x,y) ^ (150,-28) (/i = 0.99): a 40 min 
Ca II 854.2 nm time sequence started at 08:24 UT, and a 



24 min Ha sequence started at 09:15 UT. The following 
wavelength sequences were used. For the Ca II 854.2 nm 
data set we used a 29 line position sequence, ranging from 
—200 pm to +200 pm from line center (±68 km s~^), with 
10 pm steps in the core and 20 pm steps in the wings. 
12 exposures per wavelength position were recorded, re- 
sulting in a time of 11 s to complete a full wavelength 
scan. For the Ra data set we used a 25 line position 
sequence, with 10 pm steps ranging from —160 pm to 
+80 pm from line center (—73 to +37 km s~^ Doppler 
shift). At each wavelength position, 8 exposures were 
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recorded. The time to complete a full wavelength scan 
was 6.7 s. 

The image quality of the time s equences benefited from 
the SST adaptive optics system ([ Scharmer et al.| | 2QQ3b[ ) 
and the image restoration technique Mul ti-Obje ct Multi- 
Frame Blind Deconvolution (MOMFBD, |van Noort et al.^ 
2005 ) . All images from all wavelength positions m a scan 



were divided in overlapping 64x64 pixel subfields and 
processed as a single MOMFBD restoration. In such 
a restoration, the wide-band channel served as anchor 
for the narrow-band CRISP exposures to ensure precise 
alignme nt between the restored narrow-band images. We 
refer to van Noort & Rouppe van der Voort (2008) for 
more details on the MOMFBD processing strategies on 
similar extensive mult i- wavelength scans. 

After MOMFBD reconstruction of the individual line 
scans, the images from the different time steps were com- 
bined to form time series. The images were de-rotated 
to account for diurnal field rotations, and aligned and 
de-stretched to account for translation and warping due 
to seeing effects. (Local) offsets were determined on 
the wide band images and applied to the corresponding 
CRISP images. 

The Ca II 854.2 nm dataset is of excellent quality, with 
most images close to the diffraction limit of 0.23 arcsec. 
The Ha dataset is not as good, as the seeing was de- 
grading in quality and there were more frequent short 
moments of blurring, but still at least half of the images 
are close to the diffraction limit of 0.17 arcsec. 

We obtained a context Stokes V magnetogram from a 
CRISP scan of the Fe 630.25 nm line recorded at 09:39 
UT, after the Ca II 854. 2 nm and Ra sequenc es. We 
followed the procedures of "Scharmer et al. (2008) except 
that we did not use MOMFBD restoration on the polari- 
metric observations. In the lower left panel of Fig. |2j a 
blue wing (—4.8 pm) Stokes V magnetogram is shown. 
The magnetogram is the result of flatfielding and adding 
of in total 32 exposures - 8 exposures per liquid crystal 
state. We estimate the noise level to be 2.6 • 10~^ for 
Stokes V relative to Icontinuum' spatial reso- 

lution better than 0^5. The magnetogram was aligned 
to the time sequences using the wide-band channels. As 
the magnetogram was observed after the time sequences, 
no detailed information on the individual magnetic ele- 
ments can be obtained but it serves the purpose of giving 
general information on the distribution of the different 
magnetic polarities over the FOV. 

We searched both datasets for RBEs. In blue- wing 
images they appear as roundish or elongated dark struc- 
tures on top of the photo spherically formed background 
(jLeenaarts et al.||2006b|a|). In order to efficiently explore 
the datasets, we used widget based analysis tools pro- 
grammed in the Interactive Data Language (IDL) : CRIsp 
SPectral EXplorer (CRISPEX) and Timeslice ANAlysis 
Tool (TAN AT). These tools allow simultaneous display of 
spectra, images and xt-slices of arbitrarily shaped curves. 
This way it is easy to distinguish RBEs from intergran- 
ular lanes, dark reversed granulation and, by visual in- 
spection of their spectral signature, of Shockwaves. 

Our method to find and study RBEs in our data is dif- 
ferent from the Langangen et al. ( 2008[ ) approach. They 
identified RBEs purely on a spectroscopic basis, search- 
ing for specific spectral signatures in At-dataslices. We 
find numerous examples of short-lived, narrow, and elon- 



gated features in our data at wavelength offsets corre- 
sponding to significant Doppler velocities. These features 
emerge from regions with magnetic field concentrations 
and have a similar spectral signature in At-dataslices as 
the Langangen et al.' f2008) RBEs. 

To expand on the statistical properties of the RBEs 
we augmented the above described manual detection 
scheme with an automated algorithm. This algorithm 
finds RBEs by isolating the location and time of long 
thin features at high blueshifts (respectively 30 and 60 
km s~^ for the Ca II 854.2 nm and Ha datasets). To 
reject intergranular lanes and other non-RBE associated 
features, we impose that the features not be too curved. 
This automated procedure identified 413 features in the 
Ca II 854.2 nm dataset and 608 in the Ha dataset. We 
confirmed the validity of the automated search results 
by visually inspecting blue wing images and At-dataslices 
for the fiagged RBEs. This inspection confirmed that the 
automatically detected featu res have spectral c harac ter- 
istics similar to the RBEs of Langangen et al.l (12008). 

In addition to these events, we also (manually) found a 
number of small, roundish features of h igh b lueshift that 
the algorithm failed to detect (see Sec. 3.2). 



3. RESULTS 

Figure [l] shows details of one selected RBE observed 
in Ha. The bottom panel shows part of a blue- wing 
Ha image at an offset from line center equivalent to a 
Doppler shift of —59.4 kms~^. The narrow dark streak 
in the center is identified as an RBE. At this wavelength 
position, the first sign of this RBE appeared 40 s ear- 
lier, close to the photospheric bright point at the lower 
left of this feature. During its lifetime, the streak moves 
away from that bright point and vanishes after 54 s. The 
white cross marks the spatial position where a detailed 
spectrum is showed in the middle panel, and the spec- 
tral/temporal evolution in the top panel. The detailed 
spectrum is compared to the average Ha profile: the 
dashed line represents the difference of the two profiles. 
The RBE spectrum displays significant asymmetry to- 
wards the blue wing, differing from the average profile 
even beyond a Doppler velocity of 72 kms~^. 

The top panel shows the temporal evolution of the 
spectra in the central part of an RBE in a At-dataslice. 
This neatly illustrates why these events are called "rapid 
blue excursions" as they suddenly develop a highly asym- 
metric line profile towards the blue. At this spatial loca- 
tion, there is no spectral signature pointing to any down- 
fall phase associated with this event (i.e., features occur- 
ring on the red side of the line). 

Figure [2] shows the full field-of-view of the datasets 
of June 15. The upper- left panel shows a blue- wing 
Ha image at a Dopplershift of —45 km s~^. It shows 
bright network elements and many blueshifted chromo- 
spheric structures. At this Dopplershift, one can identify 
many narrow elongated features. However, some of these 
features are not identified as RBEs as they have differ- 
ent spectral signatures (i.e., not asymmetric to the blue 
wing) and/or display different temporal behavior (i.e., 
no sudden, or rapid, appearance). 

The upper-right panel shows the simultaneous Ha line- 
core image. Green curves indicate the trajectories of the 
Ha RBEs we found in an automated fashion - a dense 
concentration was found in the region of interest (ROI) 
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Fig. 3. — Detailed cutouts of ROIl (left) and ROI2 (right) showing the time-evolution of Ho; RBEs in different line positions as indicated 
above the top row of panels. A negative velocity means a blueshift. Time increases downward, with 6.7 s between successive images in a 
column. The RBEs appear as thin dark streaks moving through the field of view in the leftmost column at a doppler shift of —73.1 km s"-*^. 
The arrows indicate the RBEs when they are the most visible. They move away from the network, with their trajectory roughly aligned 
with the fibrils visible in the line-core images, but no one-to-one correspondence between fibrilar structures at the other positions in the 
line. Time-lapse animations of this figure are provided electronically. 

marked '1'. The RBEs found in R0I3 appear to have 
different behavior than the others, with movements per- 
pendicular to the magnetic field structure as outlined by 
the fibrils in the line core image. 

The lower-left panel shows the context magnetogram. 
It shows two regions of unipolar network in the lower 
left quadrant, and a small bipolar region at (x^y) = 
(35, 35) arcsec. The right side of the field shows scattered 
patches of unipolar field with opposite polarity from the 
two network regions. The lower-right panel shows the 
Ca II 854.2 nm line-core, with the Ca II 854.2 nm RBE 



trajectories overplotted in yellow. 

The RBEs are predominantly found around the two 
patches of unipolar network in the lower-right quad- 
rant of the FOV with a minority of the RBEs scattered 
throughout the FOV. The Ha-core image shows much 
more fibrilar structure than the Ca II 854.2 nm one. The 
latter shows an internetwork acoustic shock pattern (see 
e.g., Vecchio et al.|2QQ9j ) except very close to the network. 
This indicates that Ha has a much higher chromospheric 
opacity, and is likely to show structures at larger height 
above the solar surface. 
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Fig. 4. — Timeslice along the trajectory of tlie RBEs in ROIl (left-liand panels) and ROI2 (right-hand panels). They appear as diagonal 
streaks in the —73.1 km panels (they move from left to right). The —32 km s"-*^ panels show a darkening just before the RBE, extending 
less far away from the network than the RBE itself. The black dashed lines indicate the timespan displayed in Fig. [s] 



-0.15 



-0.10 



AA [nm] 
-0.05 



0.00 



0.05 



1500 



^ 1000 



500 




-60 



-40 -20 
Av [km/s] 



20 



Fig. 5. — Example line profiles of the Ho; RBEs. Black: aver- 
age spectrum over the field of view, with diamonds indicating the 
different wavelength positions of the filter. Red: profiles of RBEs; 
blue: profiles of "black beads". 

3.1. Rapid blue excursions in Ha 

Figure |3] shows blowups of ROIl and R0I2, following 
the detailed evolution of two RBEs whose trajectories 
are colored red and blue in the top right panel of Fig. [2] 
They appear as thin elongated dark streaks in the Av = 
—73.1 km s~^ column. 

We strongly encourage the reader to view the time- 
lapse animation of this Figure that is available on-line. 
Presenting the temporal evolution of these narrow, short- 
lived features on paper as a series of images cannot rival 
the visual impression of an animation. 

In the left-hand panels the RBE from ROIl is shown. 
It appears in the lower-right corner of frame 2 (numbers 
in the lower left corner). It then increases in visibility, 
becoming darker and longer while moving diagonally to- 
wards the upper-left corner of the ROI. While it traverses 
it starts to fade and is almost invisible in frame 10. The 
structure is slightly curved and follows a curved path, 
which suggests movement along a magnetic field line. 

The left-hand panels of Fig. |4] show the time evolution 
of the brightness along the trajectory of the RBE. It 
appears as a diagonal streak in the Av = —73.1 km s~^ 
panel. From this panel its apparent horizontal speed can 
be determined as 84 km s~^. 

The RBE from R0I2 is shown in the right-hand panels 
of Fig. [3] . It is present from frame 2, sticking out from 
the lower-left corner towards the center of the ROI. It 
does not appear to move until frame 4. It then moves 



towards the upper-right in frames 5-8. It fades and is 
completely gone in frame 10. The right-hand panels of 
Fig. [4] show an xt-slice of the brightness along the tra- 
jectory of the RBE. There seems to be a structure with 
similar horizontal speed in the —32 km s~^ panel (the 
black structure between x = and x = 5 arcsec). Close 
inspection of the image sequence in Fig. [3] shows that 
it is caused by a number of different fibrils and most of 
the absorption is unrelated to the RBE. The apparent 
horizontal speed of the RBE is 114 km s~^. 

The trajectories of both RBEs are roughly aligned with 
the fibrils and mottles visible in the wavelengths closer 
to the line core, but there is no one-to-one relation- 
ship between the RBEs and the dark structures in the 
Av = — 32 km s~^ indicating that the RBEs are indeed 
structures with a large line-of-sight velocity component, 
and not the damping wing of an optically thick chro- 
mospheric structure with a small velocity. There is no 
corresponding absorption in the red wing during or after 
the RBE. This rules out that our RBEs are high-velocity 
type I spicules or a related shock phenomenon. 

The two red curves in Fig.[5]show typical line profiles of 
RBEs. They are strongly asymmetric, with an extended 
absorption wing on the blue side of the core relative to 
the average profile. The largest excess absorption occurs 
between Doppler shifts of —60 and —30 km s~^. 

3.2. ''Black beads'' in Ha 

We found a large number of RBEs in R0I3. They 
are different in several respects from the RBEs discussed 
in Sec. |3.1[ Figure [6] shows an image sequence of R0I3 
with several of these special RBEs, which we dub "black 
beads". They appear as roundish darkenings, either in- 
dividually or in short strings of several beads and are 
located along a neutral line between two photospheric 
magnetic field concentrations of opposite polarity (see 
the magnetogram in Fig. [2]). Their horizontal speeds are 
relatively low (15 km s~^ j and they move parallel to the 
neutral line, but perpendicular to the large scale fibril 
orientation (see the Ha-core panel of Fig. [2|. The small 
scale magnetic field configuration close to the neutral line 
is complex, as indicated by the appearance of the line- 
core images. It is unclear how the horizontal motion of 
the black beads relates to this configuration. 

The line profiles of the black beads (see the blue curves 
in Fig. |5| are even more asymmetric than the line profiles 
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Fig. 6. — Detailed cutouts of ROI3, showing the time-evolution 
of YLa "black beads" in different line positions as indicated above 
the top row of panels. A negative velocity means a blueshift. Time 
increases downward, with 6.7 s between successive images in a col- 
umn. The white arrows indicate the black beads when they are 
most visible. A time-lapse animation of this figure is provided 
electronically. 

of the RBEs associated with the unipolar network, with 
nearly flat line profiles between Av = — 74 and Av = 
—50 km s~^. They show strong absorption all the way 
up to the most blueward position of our line scans, which 
suggests that the black beads are structures that, at least 
partly, move with an upward velocity well in excess of 
70 km s-^ 

Typically, these black beads have diameters between 
0.15 and 0.3 Mm, as measured at Av = — 74 km s~^. 

3.3. Rapid blue excursions in Ca II IR 

Figures [7) |8] and [9] show details of two RBEs in the 
Ca II 854.2 nm line, in identical format as the Ha figures. 

Figure 7\ shows the image sequences of ROM and 
R0I5. The RBEs are invisible at the largest blueshift 
column, instead they appear clearly at a blueshift of 
—30.6 km s~^. 

The left-hand panels show the appearance and disap- 
pearance of two partially overlapping RBEs. The first 



RBE appears at the arrow in the panel marked 'a'. It 
extends, reaching maximum length in in panel b, after 
which it retracts and fades. The second RBE appears in 
panel b, on top of the first one, but at a slightly different 
angle. It gains opacity in the following panels, bringing 
out its curved lower end, indicating that it is indeed a dif- 
ferent RBE than the first one. It moves upward and fades 
from view, and is gone in the last panel. The left-hand 
panels of Fig. [8] show timeslices of the intensity along the 
trajectory of the second RBE. It shows up as a diagonal 
streak in the A^' = —30.6 km s~^ panel, again without 
counterpart in the panels at other Doppler shifts. 

The right-hand panels show a slow-moving roundish 
RBE appear at the arrow in panel c. It moves slowly 
away from the network elements at the top of the ROI 
and fades from view in the last two panels. The timeslices 
in the right-hand panels of Fig. [S] show the RBE only in 
the Av = —30.6 km s~^ panel. 

Figure |9] shows typical line profiles of two RBEs. The 
profiles are asymmetric, with extra absorption in the blue 
wing, and a red wing very similar to the average profile. 
The largest extra absorption occurs around —30 km s~^, 
at significantly lower velocity than in the Ha profiles. 

3.4. Statistical properties of rapid blue excursions 

The automatic detection procedure identified 413 fea- 
tures in the Ca II 854.2 nm dataset and 608 in the Ha 
dataset. For each of the identified features we subtract 
the average spectral profile (averaged over the whole field 
of view) from the profile at each pixel along the feature. 
We use the first and second moments with respect to 
wavelength as estimates of the Doppler velocity, 'i^Doppier 
and width of the blue-shifted component, i.e.. 



^Doppler 



^^^^ 



avg 



/) dA 



/Al(^avg-/)dA 



(1) 



(2) 



where c is the velocity of light, Aq the line center wave- 
length, / the observed intensity, /avg the average inten- 
sity and Ac = Ao + Ao'^Doppier/c The integration ranges 
from the most blue- ward wavelength we observed (Amin) 
to line center, and is only executed for wavelengths for 
which /avg — / > 0, i.e., where the feature shows absorp- 
tion in the blue wing. 

We note that these estimates are necessarily rather 
crude, as it is impossible to rigorously derive the atmo- 
spheric velocity structure from the line profile. Cloud 
modelling might in some cases have yielded a more pre- 
cise estimate, but in the majority of cases the RBEs 
are overlapping a myriad of chromospheric structures, 
on which single-component cloud-model inversion is un- 
likely to yield reliable results. Based on the detailed line 
profiles we expect our estimates of '^Doppler to be on the 
low side, and most RBEs might well have components 
with a '^Doppler that is 10-20 km s~^ higher than esti- 
mated. 

The measured lengths, Doppler velocities and widths 
are shown in Fig. 10 for both datasets. The distribution 
of lengths has a sharp lower cutoff at 0.8 Mm because 
this was used as a criterion for acceptance of a feature in 
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Fig. 7. — Detailed cutouts of ROM (left) and ROI5 (right) showing the time-evolution of Ca II 854.2 nm RBEs in different line positions 
as indicated above the top row of panels. A negative velocity means a blueshift. Time increases downward, with 11.35 s between successive 
images in a column. The Ca II 854.2 nm RBEs are seen in the —30.6 km s"-*^ row, without a signature in the —61.3 km s"-*^ column. The 
left-hand panels show two RBEs. The first one appears in the panel marked 'a' from the bright point at {x,y) = (13.7,5.2) and moves 
straight up. The second one appears in the panel marked 'b' and is slightly curved. The right-hand panels show a slow-moving RBE 
appearing in the panel marked 'c'. Time-lapse animations of this figure is provided electronically. 



the automated detection algorithm. The length has been 
defined as the maximum distance along which a feature 
shows absorption in the blue wing compared to the aver- 
age spectral profile. This definition removes some of the 
arbitrariness introduced by the fact that the features are 
detected at a specific velocity. The values for the Doppler 
velocities and widths are calculated by averaging over the 
whole length of each feature. 

We find that the values for length, Doppler velocity and 
width are consistently higher for the Ha RBEs compared 
to the Ca II 854.2 nm RBEs. Ha RBEs are on average 



somewhat longer (of order 3 Mm vs. 2 Mm), have higher 
average Doppler shifts towards the blue (35 km s~^ vs. 
15 km s~^) and higher widths (13 km s~^ vs. 7 km s~^). 
We should note that there is a wide range of values for 
Doppler shifts and widths from RBE to RBE. 

In addition, the measured Doppler velocities and 
widths of the blue components vary systematically along 
the length of the RBEs. This is shown for two sample 
RBEs in Fig. [TT] We find that many RBEs show an in- 
crease of Doppler shifts and widths from their footpoint 
to the top. The footpoints (blue-coded profiles and lines) 
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Fig. 8. — Timeslice along the trajectory of the RBEs in ROM (left-hand panels) and ROI5 (right-hand panels). They appear as diagonal 
streaks in the —30.6 km s"-*^ panels (they move from left to right). There is no corresponding structure in the other line positions. The 
black dashed lines indicate the timespan shown in Fig. [t] Note the difference in the scale along the x-axis between the left-hand and 
right-hand panels. 

13 s (i.e., two time steps for Ha), up to 120 s, with peaks 
around 40 s. 

The horizontal velocity of the Ha RBEs ranges from 
km s-^ to 120 km s-^ The Ca H 854.2 nm RBEs have 
a lower maximum horizontal velocity, with most events 
having ^apparent < 40 km s"^ 

We have also measured the transverse motions of about 
35 Ha RBEs and found that most RBEs move trans- 
versely to their long axis, with transverse displacements 
of order 0.3 Mm (within a range of to 0.8 Mm) and 
transverse velocities of order 8 km s~^ (with a range from 




to 20 km s ^). This is illustrated in Fig. 15 The mea- 



-60 -40 -20 20 40 60 
Av [km/s] 

Fig. 9. — Example hue profiles of the Ca II 854.2 nm RBEs. 
Black: average spectrum over the field of view, with diamonds 
indicating the different wavelength positions of the filter. Red and 
blue: profiles of RBEs. 

typically show absorption that is at smaller Doppler 
shifts and narrower than the middle (green coded) and 
top end (red coded) of the RBEs. We see this effect for 
both Ha and Ca H 854.2 nm RBEs. This effect is vi- 
sually quite striking when we plot the measured param- 
eters for all automatically detected RBEs for a subset 
of the field of view (Fig. 12). Here we also see that the 
Ca H 854.2 nm RBEs are typically shorter and closer 
to the footpoints (photospheric magnetic field concen- 
trations) than the Ha RBEs. 

We see clear outward motions in many of the detected 
features but measuring these apparent velocities is very 
difficult to do in an automated fashion. Instead, we used 
a manual approach in which we used the extracted posi- 
tions of the features to extract xt-slices from the data 
cube (see Fig. ^ for outlines of these positions). We 
then measured the apparent velocities, ^apparent from the 
slopes of the RBEs in these xt-slices manually. This was 
difficult in most cases because of short lifetimes and un- 
clear paths in the xt-slices. We thus have ^apparent for a 
subset of only 25 RBEs for the Ha dataset and 28 for the 
Ca H 854.2 nm dataset. There is likely to be a strong 
selection effect for these measurements favoring long life- 
times and low velocities. The lifetimes of the RBEs were 
determined from the first and last occurrences in the time 
sequence. The statistics for these subsamples are shown 
The lifetimes of the RBEs range from 



sured lifetimes of these RBEs at one single height are 
between 15 and 60 s, with an average of 33 s. 

4. DISCUSSION &; CONCLUSIONS 

We have analyzed rapid blueshifted events in 
Ca H 854.2 nm and Ha data obtained with the CRISP 
instrument at the Swedish 1-m Solar Telescope. In blue 
wing time sequences, we observe these RBEs as short- 
lived, narrow streaks that move at high speed away from 
areas with magnetic field concentrations. Close to disk 
center, the Ca II 854.2 nm events we obs e rve ar e sim- 

( [20081) . We 
and Doppler velocities 



ilar to those found by [Langangen et al. 
measure similar lifetimes (45 s 

(20 km s~^). Our observations have a higher spatial and 
temporal resolution, allowing us to study their evolution 
over time and measure their apparent motion. In ad- 
dition, our automated detection algorithm allows us to 
greatly expand on the statistics of RBEs. 

The RBEs we identified in Ha have generally higher 
Doppler velocities and larger Doppler width than their 
Ca II 854.2 nm counterparts. The RBEs are longer in Ha 
than in Ca II 854.2 nm, with the Ca II 854.2 nm RBEs 
located a bit closer to the magnetic field concentrations 
that are at the root of the RBEs. For those selected 
RBEs where we could make a reliable estimate of the 
apparent velocity, we see a slight trend of higher apparent 
motion in Ha. The lifetimes of RBEs in Ca II 854.2 nm 
and Ha seem to be similar. 

suggested 



Langangen 



et al. 



(2008) 



that the 



Ca II 854.2 nm RBEs are linked to the Hinode 
Ca II H type II sp icules observed at the limb (De 



in Figs. 13 and 14 



Pontieu et al. 2007b[ ) . This suggestion was based on 
the similarity in lifetimes, spatial extent, location near 
the network and sudden disappearance or fading. In 
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Fig. 10. — Histograms of RBE properties. Length (top panel), 
measured Doppler velocity at mid-point of extracted feature (mid- 
dle panel) and width of the blue component (bottom panel) for Ho; 
(solid line) and Ca II 854.2 nm (dashed line). 

addition, the fact that RBEs exclusively show blueshifts 
corresponds well with the fact that type II spicules only 
show upward flows. 

Our observations make the connection between RBEs 
observed on disk and type II spicules at the limb much 
stronger. First of all we find the same similarities in 
lifetimes, location and temporal evolution. In time se- 
quences at high Doppler displacement from the line cen- 
ter, we also see the RBEs moving away from the network 
as narrow streaks - dynamical behavior that agrees well 
with the upward moving type II spicules. We find that 
most RBEs undergo a significant sideways motion dur- 
ing their lifet ime. This is similar to what |De Pontieu 
et al. (|2007c|) observed in type II spicules. The ampli- 
tude of the transverse velocities of RBEs (~ 8 km s~^) 
is slightly lower than that of type II spicules (of order 
12 km s~^). Our observations show a clearly separated 



blueshifted component in our line profiles which for the 
Ha line occurs at much higher velocities (of order 40- 
50 km s~^) than in the Ca II 854.2 nm line. Such high 
velocities compare well with those reported for type II 
spicules (40-100 km s~^), especially when taking into 
account the inevitable projection effects that reduce the 
line of sight velocity when observing at disk center. 

The differences in RBE properties between 
Ca II 854.2 nm and Ha can be explained by the 
difference in chromospheric opacity of both lines, with 
Ha sampling higher layers than Ca II 854.2 nm. Our 
finding of higher velocities (both apparent and in 
Doppler) in Ha as compared to Ca II 854.2 nm is 
highly compatible with the idea of RBEs being disk 
counterparts of type II spicules: the higher opacity in 
Ha allows to sample higher layers where lower density 
plasma is propelled to higher velocities. The fact that 
on average the Ca II 854.2 nm RBEs occur closer to the 
footpoints and show lower velocities than Ha RBEs fits 
in well with the Hinode/SOT limb observations which 
have shown evi dence of acceleration of plasma along 
type II spicules (De Pontieu et al. 2007b), with higher 
velocities at the top of the spicules! In addition, such 
acceleration is directly seen in our individual RBEs, 
with increasing Doppler shifts towards the top end of 
the RBEs in both Ha and Ca II 854.2 nm. 

In addition to the opacity and acceleration argument, 
the difference in velocities in Ca II 854.2 nm an d Ha 
could also par t ly be explained by the idea that |Lan- 



gangen et al. (2008) proposed. They performed nu- 
merical experiments to explore a scenario of reconnec- 
tion jets driven by constant energy release at different 
atmospheric heights. In this scenario, reconnect ion at 
lower heights (i.e., higher density) gives rise to lower ve- 
locities - this could be another reason why disk center 
Ca II 854.2 nm RBEs have lower Doppler velocity than 
the off-limb type II spicules. In Ha however, higher at- 
mospheric layers can be probed where reconnection at 
lower density would lead to higher velocities. 

More generally, we know that the limb observations in- 
troduce an observational bias towards those events that 
rocket high above the chromospheric fibrilar junk at 
lower heights. The large line of sight obscuration close 
to the limb may well, to a large extent, hide the lower 
velocity type II events at lower inclination with the sur- 
face. The highest velocities (~ 100 km s~^) of the type II 
spicules may thus be missed on the disk because of this 
selection effect, in addition to the low opacity of the low 
density tops of the spicules. The latter effect may also 
be the reason for the slightly lower transverse velocities 
we find in the Ha RBEs compared to those of the type II 
spicules at the limb. Generally, the RBEs likely sample 
slightly lower regions in the Ca II spicules where both 
longitudinal and transverse velocities are also a bit lower. 

These and other selection effects imply that the his- 
tograms of Figs. [13] and [14] should not be interpreted as 
"the" properties of RBEs. After all, it is quite possible 
that our selection criteria might exclude events that have 
too little absorption or absorption at different Doppler 
velocities. In addition, the velocities we measure here are 
only rough estimates of the real mass motion. For exam- 
ple, the apparent velocities can sometimes be a bit more 
difficult to measure and interpret. This is because visual 
inspection of the movies at various wavelengths suggests 
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Fig. 11. — Properties as function of position of an RBE in the Ca II 854.2 nm dataset (top panels) and the Ho; dataset (bottom panels). 
Color coded measured Doppler velocity (left panels), width (middle panels) and mean spectra over three portions of the RBE (right panels: 
closest to foot point in blue, middle part as dashed green, part furthest away in red and mean spectrum in black. Upper curves show the 
spectral profile, the lower curves show the subtraction of the average spectrum and the spectral profile). The extracted RBE is shown with 
a thin colored line in the left and middle panels with blue for the third of the length closest to the foot-point, green for the middle part and 
red for the part furthest away from the foot-point. The measured parameter is shown displaced to the left of the RBE. The background 
image is at line-center (middle panels), at a blue position of —36 km s"-*^ (bottom left panel) and a Dopplergram at 30 km s~-^ (top left 
panel). At both the bottom and top end of the RBEs the amount of blueward absorption drops significantly. When this happens we set 
the Doppler velocity and width to zero. This is why the black color coding at both ends of the RBEs is not indicative of lower velocities 
or widths, but rather signifies the spatial extent of the blue absorbing feature. 



that measurements of such projected velocities can some- 
times depend considerably on the wavelength studied. 
However, the presence of a separate blue-shifted absorb- 
ing component certainly indicates strong mass flows of 
at least 30-50 km s~^ in these features. 

The events we measure thus have a wide variety of 
properties (lifetimes, velocities, temporal evolution, loca- 
tion, acceleration, etc.) that are highly similar to those of 
the type II spicules observed in Ca II H on the limb. We 
are therefore confident that we have found the on-disk 
counterparts of type II spicules. 

This finding is further strengthened when we consider 
the occurrence rate of RBEs and compare it to the oc- 
currence rate of type II spicules at the limb. The latter 
can be deduced from a visual analysis of Hinode/SOT 



data such as that studied by |De Pontieu et aT] ( |2007b[ ). 
There we find between 1.5 and 3 Ca II H type II spicules 
per linear arcsec along the limb in a coronal hole. To 



compare this to our disk center observations of a coronal 
hole, we use our automated algorithm on the Ha images 
at —35 km s~^ (the peak of the RBE Doppler shift dis- 
tribution. Fig. 10) and find on average about 40 RBEs 
per image. TheHeld of view is about TO^'by 70" . Obser- 
vations at the limb show a large amount of line of sight 
superposition since spicules are quite tall (up to 10,000 
km): we can still discern spicules (of 5,000 km height) 
at the limb even when they occur up to 115 arcsec in 
front or behind the limb. This means that a typical limb 
observation samples spicules over a spatial range of 230 
arcseconds along the line of sight. This implies that our 
40 RBEs per image would translate to 40 x 230/70 = 131 
RBEs per 70 arcsec along the limb, i.e., 1.9 RBEs per lin- 
ear arcsec. This compares very favorably with the 1.5 to 
3 type II spicules we observe on average in Hinode/SOT 
limb data. 

The discovery of the disk counterpart of type II spicules 
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Fig. 12. — Properties as function of position of all detected RBEs in the Ca II 854.2 nm dataset (left panels) and the Ho; dataset (right 



panels). Color coded measured Doppler velocity (top panels) and width (bottom panels) with a background image from Av- 
(top left), A-u^— 46 km s"-*^ (top right) and line center (bottom panels). 



-27 km s-i 



provides exciting new avenues towards resolving several 
of the major unresolved issues of these features. By pro- 
viding a top view that is unhampered by the enormous 
line of sight superposition at the limb, RBEs will allow 
for a much improved study of the formation mechanism 
that drives these highly dynamic jets. Our preliminary 
finding that many of the RBEs show an increase in ve- 
locity and width of the blueshifted absorbing component 
as one travels along the jet provides a much needed con- 
straint on theoretical models. For example, such an in- 
crease in velocity along the jet may well be compatible 
with (instantaneous) velocity profiles alo ng reconnection 



jets i n 2D MHD simulations (see, e.g. [Heggland et al 
|2009|. The increase in width of the absorbing compo- 
nent along the feature, and the sudden disappearance 
or fading of RBEs at the end of their lifetime are highly 
compatible with a scenario that includes ongoing heating 
of the plasma along the jet. Such heating provides strict 



constraints on the acceleration and heating process that 
drives these jets. Detailed comparisons of the statistical 
properties of RBEs with advanced numerical models will 
be required to ascertain whether reconnection or other 
formation mechanisms (e.g., energy deposition at chro- 
mospheric heights from electron beams) can produce the 
details of the observed velocity and width profiles. 

Now that we have established the disk counterpart of 
type II spicules, we will also be able to study in detail 
the relationship between these jets and the magnetic field 
concentrations with which they are associated. Such a re- 
lationship is bound to be quite complex, since we observe 
RBEs not only in association with mostly unipolar net- 
work in coronal holes, but also in quiet Sun, and highly 
unipolar plage regions. If reconnection plays a role in 
the formation of these jets, it is likely that component 
reconnection (e.g., at tangential discontinuities) plays a 
significant role (see e.g., |Parker||1994p . By studying the 
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Fig. 13. — Histogram of RBE properties for the Ho; subsample. 
Lifetime (top panel), measured Doppler velocity at mid-point of 
extracted feature (middle panel) and apparent velocity along the 
length of the RBE (bottom panel). 

underlying magnetic field concentrations and its interac- 
tion with the flowfield, as well as the spectral line profiles 
at the footpoints of the RBEs, we should be able to con- 
strain the various suggested driving mechanisms. 

This is exemplified by the discovery of the black beads. 
These roundish features with strong absorption in the 
blue wing have many properties of RBEs with some sig- 
nificant differences. They appear along a neutral line be- 
tween two regions of opposite polarity, and not around 
unipolar network like the other RBEs. Their shape and 
high velocities suggests that these beads may be jets 
in the direction of the line of sight. This can help ex- 
plain the extremely high Doppler velocities in excess of 
—74 km s~^ in the beads. It is tempting to speculate that 
the different magnetic topology also plays a role in the 
larger velocities. Perhaps the reconnection events release 



Fig. 14. — Histogram of RBE properties for the Ca II 854.2 nm 
subsample. Lifetime (top panel), measured Doppler velocity at 
mid-point of extracted feature (middle panel) and apparent veloc- 
ity along the length of the RBE (bottom panel). 

more energy because of opposite-polarity cancellation in- 
stead of reconnection at a tangential discontinuity? A 
preliminary study of the line profiles of the black beads 
indicates that some profiles are broadened on the red side 
of the spectral line as well. This could be indicative of 
heating at the footpoint of the jet. We will study these 
events and more generally the spectral line profiles at the 
footpoints of RBEs in a more extensive follow-up paper. 

The presence of RBEs on the solar disk also provides an 
exciting way of establishing the connection of these jets 
to the recently discovered upfiows that are seen as faint 
asymmetric profiles of EUV and UV emission lines which 
are formed at transition region and coronal temperatures 
(De Pontieu et al. 2009). These upfiows in the transition 
region and corona have velocities of order 50-100 km s~^ 
and have been connected to so-called straws: jet-like fea- 
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Fig. 15. — Histogram of transverse motions of RBEs. Transverse 
displacement (top panel) and transverse velocity (bottom panel) of 
35 Ha RBEs. Most RBEs undergo significant motion perpendicular 
to their long axis. 

tures observed as highly dynamic elongated brightenings 



De Pontieu et al. (|2009l). 
withHinode and CRISP 



on the disk in the Hinode/SOT Ca II H pas sband. Since 
the observations of De Pontieu et al. (2009) lacked chro- 
mospheric velocity information, significant uncertainties 
remain about the heating of type II spicules to coronal 
temperatures. The velocities, thermal evolution and oc- 
currence rate of the RBEs we re port on here are certainl y 
compatible with the scenario of ] 
Follow-up observational studies 
will be invaluable to fully determine the thermal evolu- 
tion of these spicules and the role they play in providing 
the corona with hot plasma. 
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